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Abstract 

The latest measurements of rare b ^ s decays in the LHCb experiment have led to results in tension 
with the predictions of the standard model (SM), including a tentative indication of the violation of 
lepton flavor universality. Assuming that this situation will persist because of new physics, we explore 
some of the potential consequences in the context of the SM extended with the seesaw mechanism in¬ 
volving right-handed neutrinos plus effective dimension-six lepton-quark operators under the framework 
of minimal flavor violation. We focus on a couple of such operators which can accommodate the LHCb 
anomalies and conform to the minimal flavor violation hypothesis in both their lepton and quark parts. 
We examine specifically the lepton-flavor-violating decays B —)■ Bg 4>££', B —)■ {'K,p)££', and 

Bd,s U', as well as Kl —^ ey and K —)■ irepL, induced by such operators. The estimated branching 
fractions of some of these decay modes with yr in the final states are allowed by the pertinent experi¬ 
mental constraints to reach a few times 10“^ if other operators do not yield competitive effects. We also 
look at the implications for B —)■ and K —)• 'kvv, finding that their rates can be a few times 

larger than their SM values. These results are testable in future experiments. 
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I. INTRODUCTION 


The recently acquired data on a number of observables in 6 —)■ s decays have revealed some 
intriguing tensions with the expectations of the standard model (SM). Specihcally, last year the 
LHCb Collaboration determined the ratio of branching fractions of the decays —)■ 
and -B+ —?■ to be Rj^ = 0.745lQ;Q74(stat) ± 0.036(syst) for the dilepton invariant mass 

squared range of 1-6 GeV^. This result diverges from the lepton universality in the SM by 2.6cr. 
Moreover, earlier LHCb reported a local discrepancy at the 3.7a level from the SM prediction 
for one of the angular observables in the decay —)■ K*^. This disagreement has persisted 

after an updated analysis was done using the full LHCb Run I dataset . In addition, the latest 
measurements by LHCb 0 of the branching fractions of several rare 6 —?■ s decays favor values 
less than those estimated in the SM. 

Although the statistical signihcance of these anomalies is still too low for a dehnite conclusion, 
they may be hinting at the presence of physics beyond the SM. Subsequent model-independent 
theoretical works have in fact shown that new physics (NP) could resolve the tensions @-0.^ In 
particular, NP contributing via dimension-six operators of the form 

£sm+np = + H.c. (1) 


can produce one of the best hts to the b ^ s data if the Wilson coefficients Cf = 
contain NP effects mainly in which satisfy the condition 


—0.5. In this 


Lagrangian, and Gp denote the usual hne structure and Fermi constants, Ids,t6 are elements 
of the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix, = (1 — 75)/2, and at the rrib 
scale ~ ~ 4.2 universally for all charged leptons. In contrast, dimension-six lepton- 

quark operators with tensor structures are excluded by the measurements and (pseudo)scalar 
operators can explain them only with much hne-tuning |P], while the counterpart of £sm+np 
with right-handed quark chirality leads to much poorer hts to the data [0, . 

In view of the lepton nonuniversal nature of G^’^^ and the size of Gf’^^ relative to Gf^, if the 
tentative indications of NP are substantiated by upcoming experiments, one generally expects 
that there can be 5 —?■ s transitions which violate lepton-havor symmetry and have rates within 
reach of searches in the near future |]T|.^ Such a possibility for lepton havor violation (LFV), 


and other LFV phenomena that could have connections to the interactions of concern, have been 
examined further in the literature in the contexts of various NP scenarios 


13, 14 


In this paper, we also take these anomalies in 5 —s data to be due to NP and explore some 
of the potential consequences for a variety of rare meson decays with LFV. To do so, we adopt 
the framework of so-called minimal flavor violation (MFV), which is based on the hypothesis that 
Yukawa couplings are the only sources for the breaking of flavor and CP symmetries ||^ pTl , 
as flavor-dependent quark interactions beyond the SM are empirically ruled out if they cause 
substantial flavor-changing neutral currents. Although the implementation of the MFV principle 


^ The tensions may also be alleviated by including the effects of charm-anticharm resonances [|^. 

^ The general arguments of regarding lepton-flavor violation do not hold for certain models, such as those 
studied in |l^ . 
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for quarks is straightforward, there is no unique way to extend it to the lepton sector, as the 
SM alone does not accommodate LFV and it is still unknown whether neutrinos are Dirac or 
Majorana particles. Since there is now compelling evidence for neutrino masses and mixing [IT? 


it is interesting to formulate leptonic MFV by incorporating ingredients beyond the SM that can 
account for this observation Hi- Thus, here we consider the SM expanded with the addition 
of three heavy right-handed neutrinos as well as effective dimension-six quark-lepton operators 
with MFV built-in. The heavy neutrinos participate in the usual seesaw mechanism to endow 
light neutrinos with Majorana masses. We will focus on a couple of such operators which can 
bring about the NP contributions mentioned earlier and satisfy the MFV criterion in both the 
quark and lepton sectors.^ We will examine how these operators may contribute to a number of 
rare 6 —?■ s and b ^ d decays with LFV as well as the rare kaon decays Kl —t ep, K —)■ 7rep, 
and K —)■ ttvv. 


In the next section, after briefly reviewing the MFV framework, we introduce the NP operators 
of interest. In Section 03 we write down the decay amplitudes and proceed with our numerical 
analysis. We provide our conclusions in Section m Additional information and lengthy formulas 
are relegated to an appendix. 


II. OPERATORS WITH MINIMAL FLAVOR VIOLATON 

In the SM supplemented with three right-handed Majorana neutrinos, the renormalizable 
Lagrangian for fermion masses can be written as 


- \ {M, 


^u)kl^l,R^l,R + H.C. , 


( 2 ) 


where summation over k,l = 1,2,3 is implicit, Yu^d,u,e are matrices for the Yukawa couplings, 
Qk,L {Lk^L) denote left-handed quark (lepton) doublets, Ui^r and Di r ^ and Ei r^ represent 
right-handed up- and down-type quarks (neutrinos and charged leptons), respectively, H stands 
for the Higgs doublet, H = iT 2 H* with T 2 being the second Pauli matrix, Mj, is a matrix 
for the Majorana masses of Ui and ^ = {i'k,RY, the superscript referring to charge con¬ 
jugation. With the nonzero elements of Mi, chosen to be much greater than those of vY,y/\/2, 
the seesaw mechanism becomes operational [^, which leads to the light neutrinos’ mass ma¬ 
trix nil, = ~('i^^/2) = ^PMNs^y where v ~ 246 GeV is the Higgs’s vacuum ex¬ 

pectation value, t/pMNs denotes the Pontecorvo-Maki-Nakagawa-Sakata (PMNS |^) matrix, and 
= diag(m^, m 2 , rrtg) contains the light neutrinos’ eigenmasses m ]^2 3 - This allows one to pick 
the form 


^ i\/2 

V 


( 3 ) 


where O is in general a complex matrix satisfying 00"^ = 1, the right-hand side being a 3x3 
unit matrix. 


^ Various scenarios of leptonic MFV have been discussed in the literature M-M . 
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We assume that the right-handed neutrinos are degenerate in mass, 


= M. diag(l, 1 , 1 ) . 


( 4 ) 


The MFV hypothesis |1l6| , p! 8 | then implies that is formally invariant under the global flavor 
group = GgXG^, where G^ = SU(3)q xSU(3)f/xSU(3 )d and Gi = SU(3)l x0(3)^xSU(3 )e. 
This entails that Qk,Li Uk^R, Dk^R, Lk^L, ^k,R, and Ek^R belong to the fundamental representations 
of their respective flavor groups, 


Ql ^qQl^ 

Lr ^ VrL, 


Ur ^ 

^R U^u^Ri 


Ur -t V^Dj^, 
Ur VrUr , 


( 5 ) 


where Vq^u,d,l,e G SU( 3 )q,; 7 ^£),l^e and G 0(3)[/ is an orthogonal real matrix |]^ ^ p!9 |. 
Furthermore, under the Yukawa couplings transform in the spurion sense according to 


n ^ KsuF , 


n ^ VgY.vi , 




( 6 ) 


Taking advantage of the symmetry under Qj-, we can work in the basis where 


Y, = 


V 2 


diag( 






W = 


V 2 


diag( 


mg, m^, m. 


-) 


( 7 ) 


and the fermion helds Uk, Dk, i'k,L, i^k,R, and Ek refer to the mass eigenstates. More explicitly, 
{Ui, 1 / 2 , 113 ) = {u,c,t), {Di, D 2 , D 3 ) = {d,s,b), and (i^i, F' 2 , F' 3 ) = {e,^,T). We can then express 
Qk,L, Lk,L, and Y„ in relation to the CKM matrix I^km and 17 pmns as 


Qk,L — 


[Vckm] 

D, 


'k,L 


Uk,L — 


{UpMNs) kl ^l,L 
Uk,L 


Y.= 


V 2 


T^ckm diag(m„, m^, mj . ( 8 ) 


To put together effective Lagrangians beyond the SM with MFV built-in, one inserts products 
of the Yukawa matrices among the relevant helds to construct ^^-invariant operators that are 
singlet under the SM gauge group [|l 6 i Of interest here are the matrix products 


A, = YJ^^ = V^^^diag{yl,yl,y^)Vc 


A, = Y.Y) = 


If 

9 '-^PMNS ’ 


B. 


B. 


= = diag(|/^,2/g^|/3 , (9) 

= = diag(|/2,?/J,?/2) , (10) 


where y^ = \/2m^/v. Since the biggest eigenvalues of and Bg are, respectively, y^^l and 
y"^ ^ 3 X 10“^ at the hF-boson mass scale, for our purposes we can devise objects containing at 
most two powers of Ag and neglect contributions with Bg, as higher powers of Ag can be related to 
lower ones by means of the Cayley-Hamilton identity As for A^, to maximize the NP effects 
we assume that the right-handed neutrinos’ mass Ai is sufficiently large to make the maximum 
eigenvalue of A^ equal to unity, which fulhlls the perturbativity requirement Hence, as 

in the quark sector, we will keep terms up to order A^ and drop those with B^, whose elements 


are at most y^ 


M X 10 ^. Accordingly, the pertinent building blocks are 
= CqI + Cl Ag -I- C 2 Ag , \ = CqI + Cl A^ -I- .^2 A^ , 


( 11 ) 
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where in our mo del-independent approach (q 12 ^012 parameters expected to be 

at most of 0(1). Hence one or more of them may be suppressed or vanish, depending on the 
underlying theory. Since these parameters have negligible imaginary components pi] , |25| , we can 
make the approximations Aj = and Aj = A^. 

It follows that the ^j-invariant dimension-six operators which are SM gauge singlet and of 
the type that can readily give rise to the NP terms in Eq. (j^) are 


~ 1 


02 = 


QLlr,\2^aQ 


L '^a ) 


( 12 ) 


where the objects Agj and A^j are, respectively, of the same form as A^ and A^ in Eq. (PTj), but 
have their own independent coefficients and and a = 1,2, 3 is implicitly summed over. 
These generalize operators that were previously introduced under the assumption of MFV only 
in the quark [j^ or lepton [j^ part. The MFV effective Lagrangian of interest is then 




(13) 


where the mass scale A characterizes the heavy NP underlying these interactions. In general, one 
could also consider dimension-six (pseudo)scalar or tensor operators, or operators with a right- 
handed quark current, that fulhll the MFV requirements and contribute to 6 —?> si£' transitions. 
However, since as discussed in Section | the current data do not favor such types of operators, 
their contributions can be neglected. 

Before going into further details, it is instructive to go over the flavor structures of and (P 2 . 
Expanding them in terms of the upper and lower components of the left-handed doublets and 
and suppressing the gamma matrices, we arrive at 


“ {^ql)mS^^l)kl{UmUn^k^l + U^U^E^Ei) + {Agi)^^{An)+ D^D^E^E,) , 
^2 = {^q 2 ) ^^{^12) kl{UmUn^k^l ~ U^U^Ej^Ej) — [Ag 2 ) A ^2) kliDm^n ^k^I ~ ^^D^Ef^E^ 

+ 2(V'J'j^M^g2)^„(A£2)A:Z D^U^Uj^Ei + 2(Ag2l4j^M)mn^^^2)fcZ U^D^Ej^Ui , (14) 

where A^j = = Qjl + Cy diag( 0 , 0 , i/^) + ( 2 ^ diag( 0 , 0 , 2 /^^) with the terms 

having been dropped, summation over k,l,m,n = 1,2,3 is implicit, and = {UpMNs)kn^n 
represents a flavor eigenstate. Hence, in the approximations we have made, £mfv does not cause 
flavor-changing transitions between up-type quarks. In contrast, flavor changes among down- 
type quarks can occur with either charged leptons or neutrinos being emitted, but according 
to Eq. (j^ the two operators contribute differently to the two types of processes, which will be 
treated in more detail below. We will especially deal with the exclusive decays of & ^ q££' and 
b —)■ quu' for q = s,d and £ 7 ^ £', as well as s —?■ d££' and s —?■ dvE. 


III. DECAY AMPLITUDES AND NUMERICAL ANALYSIS 


In the presence of £mfV) one can generalize Eq. (jlj) to 

Q^e ^qb ^F 


^hqW — 


TT 


C,,,qj^PLb£^PJ' + H.c. , 


(15) 
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( 16 ) 


where is the CKM factor, q = s,d, and 


C — A - I - r 

'^ee — ^ee *-^9 '-n' > 

having set Cf,^ = —Cf^. Accordingly, for ££' = E^Ei one can derive from £mfv 

^e,e, = 4‘A, +5®^, = d'5« + y(A,)„+#(A?)„, 

where belongs to Oj and C(ii ^2 given by 


( 17 ) 


Cr ~ 'brl A ^r2 ’ ^ — 0,1,2, — 




\/2a A2G', 


J = 1,2, (18) 


the terms with in having been dropped. It follows that Based on Eq. (pi^f), 

one can write down the corresponding Lagrangian Csmi' for s —?> d£~£'^. 

From the foregoing, we can derive the contributions of £mfv to the amplitudes for a number 
of 6 —)■ s£~£'^ transitions with £ ^ £'. Thus, with = c^^,, for B —)■ K£~£'^ we obtain 






x/2 


TT 


(Pb + 


m 


B 


m 


K 


(Fo - Fi)p^ 


h,PL^' , 


(19) 


where ^ and p = Pi + p^/ denote the four-momenta of the mesons and dilepton, respectively, 
and Eqi stand for the form factors in the hadronic matrix element which 

are described in Appendix ^ and depend on the Lorentz-invariant s = p^. For B —)■ K*£~£'^ 
we arrive at 


■^B^K*U' 


^e^sb ^F 


V 2 


TT 


A ^*^PbPk* -iEe*^ + m)e*-p {p^ +Pk>), 


i'M.e*-pp^ 


(Y'PlI', 


( 20 ) 


A = 2V/{niB + rriK*), C = Ai (tub-|- m^*) , ®) = A 2 /{rriB + mx*) , 

sH = C — D — 2^0771^* , (21) 


where V and Ao,i ,2 are the form factors for {^^*187^(1—75)611?), which are defined in Appendix 
We look at Bg —)■ (j)£~£'^ as well, which has an amplitude -Mss^ipii' analogous to M.B^k*u'- 
Bg £-£'+ we find 


-^1+ [mi' + , (22) 

where is the Bg decay constant, which is also dehned in Appendix]^. 

The MFV Lagrangian in Eq. (|^ generates lepton-flavor-violating b —)■ d£~£'^ processes with 
the same coefficient c^, in Eq. (piBf), but a different CKM factor, This makes it of interest to 
include them in our study, particularly B ti£~£'^ , B p£~£'^ ^ and Bd £~£'^. Generally, 
there could be other dimension-six MFV operators that also contribute to 6 —?■ d£~£'^ and 
therefore can enhance or reduce the impact of c^^,. Hereafter, we focus on the possibility that 
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the effects of such operators are unimportant. Under these assumptions, the rates of these 
decay channels have expressions similar in form to those for B —)■ B —)■ 

and Bg —)■ respectively, but are comparatively smaller because \\db/Kb\'^ ~ 1 / 22 . 

Before starting our numerical calculation, we need to specify our choices further. Since the Ai 
matrix in Eq. (pi]|) can be realized in many different ways, we concentrate on the least complicated 
possibility that O is a real orthogonal matrix, in which case 


— 


2M 


u TV r/1 

’^PMNS '-^PMNS 


(23) 


For UpMNSi we employ the standard parametrization o , with its elements being determined from 
the results of a recent fit to global neutrino data in Ref. @] , which depend on whether neutrino 
masses have a normal hierarchy (NH), < mg < m 3 , or an inverted one (IH), m 3 < m^ < mg. 

Since the empirical information on the absolute scale of mi 2,3 is still far from precise for 
definiteness we pick mi = 0 (m 3 = 0) in the NH (IH) case. Requiring the biggest eigenvalue of 
to be unity then implies ~ 6.1 x 10^“^ GeV. For the elements of V/km, we adopt the results 
of the latest fit performed in Ref. |27|. With these numbers, we can evaluate the branching 


fractions of the decay modes discussed earlier using the rate formulas collected in Appendix 
which also describes our choices for the relevant hadronic form factors and decay constants. 

We now attempt to attain the largest branching fractions of the b —)■ q££' decays of interest, 
for I ^ ^ under our MFV framework by scanning the space of the parameters, which enter 

the NP term, c^^,, in the Wilson coefficient according to Eqs. (pUD-(p! 8 D. This amounts to 
maximizing |c£^/|. Simultaneously we need to impose the pertinent restrictions from the existing 
6 -meson data. Thus, based on the results of an analysis of the latest b ^ s measurements 
including the LHCb anomalies [|^ Q, we require 


c„ = 0 , -0.71 < < -0.35 , (24) 

which can lead to one of the best fits to the data p. Nevertheless, it is possible to let Cgg have 
some nonvanishing value as well [0. Accordingly, we alternatively impose 


0 


< Cee < 0.3 , 


-0.65 < c 


74/4 


< 


-0.45 


(25) 


which is well within the allowed la best-ht region in the second plot of Fig. 6 in Ref. 0. 

In addition, since for a specihc pair of £ and £' ^ £ affects also s d££' processes, we 
need to take into account the available experimental bounds, which may imply complementary 
limitations on c»£,. Indeed, if like before other dimension-six MFV operators do not give rise 


to competitive effects on s —)■ d££'^ we find in kaon data that the branching-fraction limit [|W 
B{Kl —>■ e^p^)gxp = {J3{Kl e+/i“) -|- B{Kl —)■ ^ ^ 10 “^^ can translate into the 

strongest restriction on among lepton-flavor-violating meson decays. Thus, applying Eq. (| 
in the appendix, with the central values of the input parameters, we then extract 


< 0.16, 


(26) 


which we will also impose. 

In our scans of the parameter space to maximize |c££/|, we hnd that the bound in Eq. (^) 
is always reached. This implies that it can already be used to estimate the largest branching 
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fractions of various b —)■ ge/r decays and K —)■ Tre/r within our MFV scenario with Oi ^2 taken 
to be the main operators responsible. We display the results in the third column of Table | 
and compare them with the corresponding experimental limits if available."^ One observes that 
the predicted B{B —)■ is only 4 times below its measured bound and, therefore, may 

be probed in near-future searches. For the other modes, the predictions are lower than their 
experimental counterparts by more than an order of magnitude. 


For the decay channels with if = er or pr, one can entertain many different possibilities. 
For several of them, we present the results listed in Table ||, where we have separated those 
obtained under the constraint in either Eq. or Eq. (|^) , besides the requirement in Eq. 
Furthermore, to get the numbers without (within) parentheses in the table we have employed 
the central values of the neutrino mixing parameters from Ref. ^6\ associated with the normal 
(inverted) hierarchy of light neutrino masses, except for the Dirac CP-violation angle 6 which 
has a greater uncertainty than the other parameters. To reflect this uncertainty, the left and 
right numbers inside the square brackets have been computed with the minimum and maximum, 
respectively, of (5/degree = SObl^g from Ref. @ in the NH (IH) case. 


We remark that in Table || each vertical set of results for the er (pr) channels comes from 
the same maximized value of |Cer| (|c^t|)- However, for the er and pr numbers in the same 



TABLE I: Predicted upper limits on the branching fractions of exclusive meson decays involving ep in 
the final states, calculated with |ce/i| from the empirical limit B(^Kl —)• ^ ^ at 90% 


confidence level |17], under the assumption that the effects of operators Oi ,2 dominate these processes. 
For comparison, the experimental counterparts are also displayed if available. 


^ In conformity to the experimental reports |^ , the B —>■ prediction in this table is the simple average 

over the and channels, namely B[B —>■ A^*h^p^) = [B{B'^ —>• K''*''~^e^fi^)+B{B^ —>■ iy^*^°e^p^))/2, 
whereas the B —?> 7rep prediction is from B{B —>■ Tre^p^) = {B{B^ —>■ 7r+e^p^) + 2B{B^ —>■ 7r°e^p^))/2 and 
similarly for B —> pe^fi^. 
























Decay mode 

Branching fractions 

Measured 
upper limit 
at 90% CL |29| 

Predictions 

(I) 

(11) 

B+ K+e^T^ 

3.0 X 10“® 

[0.9, 2.7] ([1.5, 2.2]) X 10-® 

[1.1, 3.5] ([2.0, 2.6]) X 10-® 

B+ ^ K*+e^T^ 

- 

[1.7, 5.3] ([3.0,4.3]) X 10-® 

[2.2, 6.9] ([3.9,5.!]) x 10-® 

Bs — )• 

- 

[1.7, 5.0] ([2.8,4.!]) X 10-® 

[2.1, 6.6] ([3.7,4.8]) x 10-® 

Bs 

- 

[0.9, 2.6] ([1.5,2.!]) X 10-® 

[1.1, 3.5] ([1.9, 2.5]) X 10-8 

B~^ —>■ 7r~^e~T~^ 

2.0 X 10“® 

[0.2,0.7] ([0.4, 0.5]) X 10-9 

[0.3, 0.9] ([0.5,0.6]) X 10-9 


- 

[0.8, 2.4] ([1.4, 2.0]) X 10-9 

[1.0, 3.2] ([1.8, 2.3]) X 10-9 

B^ 

2.8 X 10-5 

[0.3,0.8] ([0.5, 0.7]) X 10-9 

[0.3,1.1] ([0.6,0.8]) X 10-9 

B~^ 

4.8 X 10-5 

[0.6,0.9] ([0.4, 0.5]) X 10-^ 

[0.8,1.4] ([0.5,0.8]) X 10-7 

B+ ^ K*+p^t^ 

- 

[1.1,1.8] ([0.7,1.0]) X 10-^ 

[1.5, 2.8] ([1.0,1.5]) X 10-7 

Bs -)• 

- 

[1.1,1.7] ([0.7,1.0]) X 10-7 

[1.5, 2.6] ([0.9,1.4]) X 10-7 

Bs p'^T^ 

- 

[0.6,0.9] ([0.4, 0.5]) X 10-7 

[0.8,1.4] ([0.5,0.8]) X 10-7 

B'^ — ^ Tr~^p^T^ 

7.2 X 10-5 

[2.8,4.6] ([1.8, 2.6]) x 10-9 

[3.9, 6.9] ([2.4,3.8]) x 10-9 

— > p'^p^T^ 

- 

[5.1,8.4] ([3.2,4.8]) x 10-9 

[7.1,13] ([4.4, 7.0]) X 10-9 

B^ — > 

2.2 X 10-5 

[1.7, 2.8] ([1.1,1.6]) X 10-9 

[2.4,4.3] ([1.5, 2.4]) x 10-9 


TABLE II: Predicted upper limits on the branching fractions of exclusive 6-meson decays involving 
(e, ij,)t in the final states, computed with the maximal determined under the imposed constraint 

set (I) Cgg = 0, —0.71 < < —0.35, and |cg^| < 0.4 or (II) 0 < Cgg < 0.3, —0.65 < < —0.45, 

and |cg^| < 0.4, as discussed in the text, under the assumption that the effects of operators Oi ^2 
dominate these lepton-flavor-violating processes. The numbers without (within) parentheses correspond 
to neutrino mixing parameters belonging to the normal (inverted) hierarchy of neutrino masses, whereas 
the left and right numbers inside square brackets reflect the minimum and maximum empirical values 
of the Dirac phase 5 in C/pmns- For comparison, the data are also displayed if available. 

column, the respective sets of values which maximize |cgi-| and |c^r| are not always the same. 
Moreover, for the er (pr) results in different columns, the sets which maximize |Cgr| (|c^t|) 
are generally also different, as we are interested in attaining the biggest branching fractions 
in the different scenarios. We note, in addition, that the results in the NH (IH) case follow 
from |Cg^|<0.9 (0.8) and |c^^|<1.9 (1.4), with |f+| ~ 0.2-0.7 (1-2) and ~ 4-17 (40-120), 
implying that A/|(CijT/t + C 2 jyt)^rj\^^^ > 1-7 (0.65) TeV if ~ 

It is evident from Table |l|that for each decay mode the results in the different cases are roughly 
of similar size and differ from each other within only factors of 4 or less. More interestingly, we 
notice that a few of the predicted branching fractions of the pr channels can be as high as a few 
times 10“^. Although they are still at least about two orders of magnitude below the existing 
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empirical limits, which are presently not many, upcoming experiments will expectedly offer ample 
opportunities to look for these decays and improve the data situation. 

Since the restraint in Eq. (^) would lessen if there were other operators having destructive 
interference with in the Ki decay amplitude, their presence would bring about a different 
set of predictions for the processes listed in Tables | and |^. Changes in the predictions would 
also occur if other operators significantly modified b —?■ dii'. Thus, our scenario in which Oi ^2 
dominate these LEV transitions will be tested when one or more of them are discovered and the 
acquired data compared to the predictions. 

It is worth commenting, in addition, that the parameter values which are responsible for 
the predictions above also translate into reductions in the rates of the lepton-ffavor-conserving 
b —>■ (s, d)T^T~ decays by up to a few tens percent with respect to their SM estimates. This 
implies that future observations of these processes with good precision will serve as important 
complementary tests on our NP scenario. 

Now, as indicated in the preceding section, £mfv in Eq. ([I^ also contributes to transitions 
with neutrinos in the final states. Specifically for b —)■ quu' decays the effective Lagrangian is 
given by 


C 


bquu' 


\/2 A-l Gt^ 

- + H.C., 

TT ' 


^uu' — ^vv' 


SM 




(27) 


where ~ 
for vu' = 


-6.4 is the SM prediction and c^, 


arises from NP. From £mfV) one then gets 


~( 1 ) ~( 2 ) 

■PkPl - ^EkEi ^EkEi 


io + tl + ?2 . 


(28) 


where enter in Eq. (pTf ), but with the opposite relative sign, and — ^^ 2 - 


Therefore, and are generally independent of each other 


Since only the contributions of with /c = / = 1, 2, 3 can interfere with the SM contribution 
to & —)■ quE, and since the neutrinos are not detected, it is straightforward to derive the ratio of 
branching fractions 


B{B /fWi/z/) 



1 + 


^kVk 

CSM 




(29) 


These channels are not yet observed, but there are experimental limits on their branching frac- 

B{B —)■ iPz/i/)g^p < A.3B{B —)■ and so we 


tions [2£l. Here the relevant bound is 


can impose 




(30) 


Their counterparts in the kaon sector, K'^ and Kl are similarly affected 

by In its presence, the effective Lagrangian for sduu' interactions is given by 


C 


\/2 OgAj Gp 


dsvu' 


'K St 


Xw + H.C. , 


^vv' — 


(31) 
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where Xt = Xsd, Sw = sin^^^w = 0.231, the SM term [|g = X + \Vus\^Re{V*^V^^)PjX^. 

consisting of 2 contributions from top- and charm-loop diagrams, respectively, with X = — 
and |K.rRe(K:Kd) = A72-A^ in the Wolfenstein parametrization, and due 

to NP. The latest predictions of the SM are B{K^ —)■ Ti^yy)sM = (9.11 ± 0.72) x 10“^^ 
and B{Kl —?■ = (3.0 ± 0.3) x 10“^^. Generalizing the expressions for these branching 

fractions to include the £mfv effects, which modify the X part, we obtain 


BiK^ 


K 


Truu = 


+ 


3A10 


BiKL vr°z/i/) 


E| [ReA, (X - + (AV2 - A'')PJ" 

+ + |A44E |c„,„,|d . 

l^k J 

+ 


Kl (ImAt)2 ^ 


3A10 




l^k 




(32) 

(33) 


where k'^ = 0.997 x 5.173 x 10-i\ X = 1.481, = 0.404, and = 2.231 x 10-^° are the 

central values from Ref. [^. In view of the data B{K^ 7r+z/iA)gxp = (17.3^Jq;5) x 10“^^ 
and [^] B{Kl vr°z/i/)exp < 2.6 x 10“®, we may then demand only 


0.7 < 


^ K+^tt+uv 


B{K^ Tl^ixp) 

B{K+ 7r+z/iA)sM 


(34) 


In addition, comparing Eqs. (p9D and (^), one infers that is equal to its —)■ 

counterpart, Tkj^^ttOuu- 


For the parameter values that yield the examples in Table |1|, if = 0 in Eqs. (|I^) and (H), 
we hnd that can be as large as 1.22 (1.15) in the NH (IH) case, whereas 


rx+^-K+vv can reach 1.15 (1.11). 


If - 




= 0 instead, the results above for the decay modes with 


charged leptons are unchanged, but now the branching fractions of the modes with neutrinos 
tend to be reduced by up to 14% with respect to their SM expectations. All of these numbers 
are well within their corresponding restrictions in Eqs. (|30D and (^). 


If both 


c^/ and 


,( 2 ) 

-w 


are nonzero, then in general and are not connected and, con¬ 


sequently, they can be maximized independently. In that case, our results in Table || are still 
the same, but for the channels with neutrinos, after scanning the parameter space subject 
to Eqs. (^OD and (^4)), we obtain a maximum value of = VK^^TrOvu that saturates its 


limit of 4.3 and vk+^-k+uv that reaches -^3.1, with r 
measurements can offer signihcant checks on these predictions. 


-7 and |^^ 2 l ~ 0. Thus future 


Finally, we mention that, in the approximations we made, /Imfv produces vanishing effects 
on the SM-dominated transitions b —)■ {u,c)ui, s —)■ uui, and c {d,s)i'i. Although its 
contribution to t ^ bui is nonzero. 


M 


MFV_ 


‘^^tb{Co2 + CuVt + C22 7)(^^2) 


kl 


A 2 




(35) 


we estimate that its branching fraction is only under 10“^ for optimistic choices of the parameter 
values, and so it is much smaller than B(t ^ bW^ —?■ bv^'^)^y[ ~ 0.1 which is consistent with 
the data |1l7| . 
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IV. CONCLUSIONS 


We have entertained the possibility that the anomalies recently detected in the measure¬ 
ments of rare & —?■ s processes are NP signals and explored some of the potential implications 
for a number of exclusive lepton-flavor-violating meson decays. Adopting the effective theory 
framework of MFV based on the SM plus 3 heavy right-handed neutrinos participating in the 
seesaw mechanism, we concentrate on a couple of dimension-six 4-fermion operators which ac¬ 
commodate the interactions that can yield one of the best hts to the 6 —)■ s data. Assuming 
that these operators conform to the MFV hypothesis in both their quark and lepton parts, we 
evaluate their effects on various meson decays that violate lepton-flavor symmetry, subject to 
additional relevant 6 -meson and kaon constraints. For simplihcation, we further assume that 
other dimension-six MFV operators do not induce competitive contributions to these transitions. 
This scenario can be tested when they are discovered and the resulting predictions confronted 
with the acquired data. 

With the preceding premises, our numerical work shows that for decays with charged leptons in 
the hnal states is the most restricted among the Wilson coefficients for £ 7 ^ its strictest 
bound being supplied by the experimental limit for —)■ e/i. Interestingly, among the other 
modes considered, the resulting prediction for B{B ^ Kejj) is the closest to its empirical bound, 
being only 4 times smaller, and therefore may be probed by forthcoming searches. Moreover, the 
predicted branching fractions of the exclusive b —)■ s^r decays can be as large as a few times 10“^ 
and hence may be examined as well in near-future experiments. For decay channels with neutrinos 
in the hnal states, B —)■ and K —)■ we hnd that they are comparatively far less 

restricted and that the impact of the MFV operators can enhance their branching fractions by 
up to a few times with respect to the SM expectations. Thus, planned measurements on these 
processes involving neutrinos will provide important checks on the enhancements. 
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Appendix A: Decay rates 


In what follows we assume that the charged leptons ^ and i' in the hnal states of the decays are 
diherent in havor. Furthermore, the decays are induced by left-handed current-current operators 
of the type in Eq. (pAf). 

The hadronic matrix elements pertinent to the decay B Kii' are then 


{K\s-f%\B) 


m 


B 


m 


^ rP F 
P ^0 


+ 


pI+pI- 


m 


B 


m 


K prj 


Fi 


(Al) 


and {K\s'y^'j^b\B) = 0, where (j^k) Pb (Pk) are the B (iP) mass and four-momentum, 
respectively, p = Pb ~ Pki and the form factors Fq ^ depend on the Lorentz-invariant s = 


12 




according to [34 


^0 ~ %o 


“ 1 “ CLq-^Z “ 1 “ Qj^cyZ CLq^Z 


^ 02 ^ 


^03^ 


= 


''10 


+ a^iZ + 0 ^ 2 ^^ ~ l(®ii ~ 2 a;i 2 )^^ 


1 - s/(mg + 0.04578 GeV)' 


z = 


— S — \/t+ — to 

— S + ■^/t^ — to 


t± = {m^+ ±mj^+y , to = (^1 - \/l - t_/t+^t+ , (A2) 


the a’s being constants. Numerically, we adopt (ctoo; ®oi) ® 02 ) ^oa) = (0.54,-1.91,1.83,-0.02) 
and (o^O) Oil) Q 12 ) = (0.43,-0.67,-1.12), which are their central values from Ref. |^. For 
the R —)■ TT form-factors, which are relevant to i? —)■ irii' and dehned analogously to Fo,i 
in Eq. (|A1|) , we make use of the parametrization choice preferred in Ref. [^^ as well as the rela¬ 
tion (7r“|d7'^5|i?“) = \/2 (7r°|(J7^6|.B°) based on isospin symmetry. 

Similarly, for B K*ii' the requisite matrix elements are 

2 R 


{k*\s^^h\B) = 


ttIb + rrij^. 


'-'rfpCTT 


£*^PbPk* 


(77*157^75615) = 2 iAQmj ^, 


e ■ p 


p'^ + iA^ (m^ + m^*) ( e 


>*77 


• P 


■P' 


iA^ e* ■ p 


m. 


pI + pI*- 


m 


B 


m 


K* 


P 


, „ , (A3) 

T V ^ J 

where p = Pb ~ Pk* the form factors V and Ao,i ,2 are functions of s = p^. For Bg —077', 
the hadronic matrix elements have expressions similar to those for B —)■ K*i£'. In numerical 
work, we adopt the B —)■ K* and 5^ —)■ 0 form-factors available in Ref. [^ from combined 
hts to light-cone-sum-rule and lattice results. For the 6 —)■ d transition B —> p77', the form 
factors are dehned analogously to those in Eq. (|A3|) , we again utilize the parametrization provided 
in Ref. and as in the 5 —)■ tt case we have (p“|d 7 ''(l —75)6|5“) = \/2 {p^\d'y^{l —75)6|5°). 

The formulas in the last two paragraphs lead us to the amplitudes in Eqs. (|^) and ( p^) , 
respectively. Subsequently, we arrive at the differential rates 

Fif2g + {ml - mlf 




^F| 




ds 




1536 Ti^m'B s 


3 a3 


(A4) 






ds 


6144 TT® ml ml, P 


2 [ (A^^s + 6 C^s + CD^) 02 + 2 Aq /g g'j ml, 
+ + 2 CD (s — ml) + D^^] 02^ 


(A5) 


where s = {Pi + PiiY, 

01 = /C(5, 

9 = F{ 

~9 = 


m^,m^, 


). 

mB,mj^,s) , 


02 = 25 ^ - (m^ + ml)s - {m\ - m^Y , 
03 = 3(m^ + ml)s — 3(m^ — m^Y ) 


)C{x, y, z) = x'^ + y'^ + z"^ — 2{xy + yz + xz) 


(A6) 


In computing the branching fractions, we integrate Eqs. ( |A^) and ( |A5|) over the whole kinematical 
ranges {m^ + m^/)^ < s < {ms — mxY ^md (m^ + m^/)^ < s < {mB — mK*Y, respectively. We 
do likewise for Bg 4>ii' and B —)■ (vr, p)ii'. 
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To examine the purely leptonic decay Bq ii', the hadronic matrix elements we need are 
( 0 |g 7 '' 6 |-Bg) = 0 and (o|^ 7 ^ 756 |i?q) = —ifs^Ps^ where is the decay constant. The amplitude 
in Eq. then follows, leading to the decay rate 




32 TT^ 


m 


B 


2 I 2 


mt 


m 


2\2 


m 


B 


(A7) 


The kaon reaction Kl —t e“/i^ proceeds from the components of ~ + K^^/\/2 both 

decaying into Analogously to Bq —)■ £i', the necessary hadronic matrix elements are 

(olsy^dlA''’) = ( 0 |d 7 ’^s|A'°) = 0 and A'°) = ( 0 |d 7 ’^ 75 s|J?°) = where is 

the kaon decay constant. Neglecting m^, with the aid of Eq. (|A7|) we then find 


T 


KL^efj. 


Ihvr^ 



(AS) 


The amplitude for Tr+e^/r’*', also containing Ce^, is similar to that for —)■ 

which arises mainly from SM interactions described by = —\/8 Gpl4*s s^^Pi^u + H.c. 

It follows that one can conveniently express the branching fraction B{K^ —)■ 7 r+e“p^) in relation 
to the well-measured —?• 7r°z/p+)exp = (3.353 ± 0.034) x 10“^, upon assuming isospin 

symmetry and neglecting rUg, without having to know the A' —)■ vr form-factors in great detail. 
Thus, since {'ii^\s'^'^d\K^) = \/2 ( 7 r°|s 7 ^M|we arrive at 

B{K+ 7r+e-p+) ~ 2 |t/T] ^ 7r°i//i+)g^p . (A9) 


Similarly, for Kl —)■ vr°e /i’*', since ( 7 r°|(i 7 '^s|A'°) = —{Ti^\s^'^d\K^) = ( 7 r°|s 7 ^M|A'’''), we get 


, |q;„c„„ ImA^jP 

B{Kl ^ TfOe^p^) = 2B{Kl ^ 7r0e-p+) ~ \ B{K+ ^ 7rVp+)g,p . (AlO) 

^ '^K+ I 


In numerical applications of these rate formulas, we employ the B^ and Bg^d {Kl and 
lifetimes from Ref. [|^ (0)’ *^e = 1/133, and Gp = 1.1664 x 10“^ GeV^. For the decay 
constants, we adopt the central values of = (190.5 ± 4.2) MeV, = (227.7 ± 4.5) MeV, 
and = (156.3 ± 0.9) MeV from Ref. ||37|] . 
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